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1 Fully packed loop models on the square lattice 

A fully packed loop (FPL) model on the square lattice is the statistical en- 
semble of all loop configurations, where loops are drawn on the bonds of the 
lattice, and each loop visits every site once [4,18]. On finite geometries, loops 
either connect external terminals on the boundary, or form closed circuits, 
see for example Figure [U In this chapter we shall be mainly concerned with 
FPL models on squares and rectangles with an alternating boundary condi- 
tion where every other boundary terminal is covered by a loop segment, see 
Figure [T] 
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Fig. 1. Fully packed loops inside a square with alternating boundary condition. 



2 Jan de Gier 



An FPL model thus describes the statistics of closely packed polygons on a 
finite geometry. Polygons may be nested, corresponding to punctures studied 
in Chapter 8. FPL models can be generalised to include weights. In particular 
we will study FPL models where a weight r is given to each straight local 
loop segment. The partition function of an FPL model on various geometries 
can be computed exactly using its relation to the solvable six-vertex lattice 
model. It is well known that the model undergoes a bulk phase transition at 
t = 2. 

We furthermore study nests of polygons connected to the boundary. In 
the case of FPL models with mirror or rotational symmetry, the probability 
distribution function of such nests is known analytically, albeit conjecturally. 
FPL models undergo another phase transition as a function of the boundary 
nest fugacity. At criticality, we derive a scaling form for the nest distribution 
function which displays an unusual non-Gaussian cubic exponential behaviour. 

The purpose of this chapter is to collect and discuss known results for FPL 
models which may be relevant to polygon models. For that reason we have 
not put an emphasis on derivations, many of which are well-documented in 
the existing literature, but rather on interpretations of results. 

1.1 Bijection with the six- vertex model, alternating-sign matrices 
and height configurations 

There is a well-known one-to-one correspondence between FPL, six-vertex 
and alternating-sign configurations [25,69]. In the six- vertex model, to each 
bond of the square lattice is associated an arrow, such that at each vertex 
there are two in- and two out-pointing arrows, see e.g. [7]. There are six local 
vertex configurations which are given in the top row of Figure [5J The six- 
vertex and FPL configurations are related in the following way. The square 
lattice is divided into two sublattices, even (A) and odd (B). For each arrow 
configuration we draw only those bonds on which the arrow points to the 
even sublattice. If we choose the vertex in the upper left corner to belong 
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Fig. 2. Bijection between six-vertex and FPL vertices. The correspondence is dif- 
ferent on the two sublattices A and B. 
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Fig. 3. An equivalent six-vertex and fully packed loop configuration. 

to the even sublattice, the six- vertex and FPL configuration in Figure [3] are 
equivalent, as can be seen from the correspondence in Figure [5] 

Alternating sign matrices (ASMs) were introduced by Mills, Robbins and 
Rumsey [51,52] and are matrices with entries in {—1, 0, 1} such that the entries 
in each column and each row add up to 1 and the non-zero entries alternate 
in sign. A well known subclass of ASMs are the permutation matrices. Let us 
also introduce the height interpretation of an ASM. Let A = (ajj)™,- =1 be an 
ASM, then define the heights hij by 

= n — i — j + 2 <Xi'j>. (1) 

i'<i, j'<j 

This rule ensures that neighbouring heights differ by one. The correspondence 
between the three objects is given in Figure 3J An example of a six vertex and 




-1 1 

Fig. 4. The six vertices and their corresponding heights and ASM entries. 



its corresponding height configuration is given in Figure[5]for the 3x3 identity 
matrix. 

1.2 Structure 

As each external terminal, or outgoing bond, is connected to another terminal, 
FPL diagrams can be naturally labeled by link patterns, or equivalently, two- 
row Young tableaux or Dyck paths. For example, the diagram in Figure [5] has 
link pattern (( () (()) () )) which is short hand for saying that 1 is connected to 
12, 2 is connected to 11, 3 to 4 etc. The information about connectivities can 
also be coded in two-row standard Young tableaux. The entries of the first 
row of the Young tableau correspond to the positions of opening parentheses 
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Fig. 5. Vertex and height interpretation corresponding to the 3x3 identity matrix. 
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Fig. 6. An FPL diagram with link pattern (( () (()) () )). 

'(' in a link pattern, and the entries of the second row to the positions of the 
closing parentheses ')'. The FPL diagram of Figure [6] carries as a label the 
standard Young tableau given in Figure [7J 
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Fig. 7. Standard Young tableau corresponding to the FPL diagram in Figure [6] 

Yet another way of coding the same information uses Dyck paths. Each en- 
try in the first row of the standard Young tableau represents an up step, while 
those in the second row represent down steps. The Dyck path corresponding 
to Figure [JJ is given in Figure [H 

In this section we collect some structural results regarding local update 
moves of FPL models. Following Wieland [77] , we define operators G„ that 
act on the height configurations as follows. They act as the identity on each 
square except on the square at where they either increase or lower the 
height by 2 if it is allowed. A change of height is allowed if neighbouring 
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Fig. 8. Dyck path corresponding to the FPL diagram in Figure[6]and the standard 
Young tableau in Figure [7] 



heights still differ by one after the change. If it is not allowed, Gij acts as the 
identity. For future convenience we also define the operators 



Go 



n g < 



n 



(2) 



where Go and G\ denote the even and odd sublattice of the square lattice 
respectively. 

Starting from an initial height configuration, such as the one in Figure [5l 
the operators Gij generate all height configurations. Put in other words, if we 
denote the height configuration corresponding to the unit matrix by Z\, all 
other allowed height configurations correspond to a word in the operators Gij 
acting on Z\. 

On a plaquette of an FPL configuration, the involution G acts as 



G 



(3) 



while on other types of plaquettes G acts as the identity. Wieland [77] observed 
that the operator Go o Gi "gyrates" a link pattern and that the number of 
FPL configurations is an invariant under gyration. 

We define two other operations on the FPL diagrams, Uij and Oy, that 
leave the link pattern invariant but that generate all diagrams belonging to a 
fixed link pattern. The operator U acts on two plaquettes, either horizontally 
or vertically. Where it acts non-trivially it is given by, 



U : 



(4) 



The operator O acts on three plaquettes, either horizontally or vertically. 
Where it acts non-trivially, it is given by, 



O 



□ D 
□ 



(•5) 



n 
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It is easy to see that both U and O leave the link pattern external to the 
plaquettes on which they act invariant. It is also not difficult to see that on a 
horizontal strip of arbitrary length, such that only the leftmost and rightmost 
edge are connected to the outside world, the operators U and O generate all 
possible FPL diagrams leaving the link pattern invariant. A similar argument 
holds for vertical strips. This proves that acting with U and O on an FPL di- 
agram with given link pattern, one generates all FPL diagrams corresponding 
to that link pattern, and no more. 



Fig. 9. An isolated row inside an FPL configuration: only the leftmost and rightmost 
edge are connected to the rest of the FPL configuration. The operators U and O 
generate all possible configurations within the row. 



2 Partition function 

To each local FPL vertex we assign a weight Wi and dehne the statistical 
mechanical partition function Z n as the sum over all FPL configurations of 
the product of the vertex weights, 

z n = yi w 

configurations i=l 

where k t is the number of vertices of type i. We will only consider the case 
where the weights on the two sublattices are the same, i.e. in the six-vertex 
representation the weights are invariant under arrow reversal. Using standard 
six-vertex notation we write Wi — W2 = a, = W4 = b and W5 — wq = c, see 
Figure Q3B 



a a b b c c 

Fig. 10. Weights of the six local FPL vertices. 
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It is convenient to parametrise a, b and c in the following way, 

a = sin(7 — it), 6 = sin(7 + u), c = sin(27), (7) 
and to introduce the w-independent quantity r by 

t2 = c 2 - (a - bf = 2 _ 4cQg2 

where Z\ is the standard notation for the anisotropy parameter of the six- 
vertex model defined by 

a 2 + b 2 - c 2 , , 

^ = ^^ = -co S (2 7 ). (9) 

When a — b, t — c/a gives a weight to straight loop segments. It is therefore 
expected that for some critical value of r there is an ordering transition in the 
FPL model from a disorder phase to a phase where the vertex with weight 
c dominates and the polygons are elongated. We will see below that this 
transition takes place at r = 2. For a > b + c or b > a + c there is another 
ordering transition at r — where the vertices with weight a or b, respectively, 
dominate. 

The partition function Z n can be computed exactly for finite n applying 
methods of solvable lattice models to the six- vertex model with domain wall 
boundary conditions. This was first done by Korepin and Izergin [36, 37, 39] 
who derived the following determinant expression for Z n , 

2 

(sin( 7 + M) sin(7-M)) n 
Ai = ; —n cr«, (1U) 



where a n is the Hankel determinant 



det -j-rrA , (11) 



Kti<t 



and 



<p(U) = r. (12) 

sm(7 -(- u) sm(7 — u) 

Using the height representation |T]) it is possible to introduce elliptic 
weights, rather than the trigonometric weights ((7|). The partition function 
in that case has been computed by Rozengren [70]. 

2.1 Another form of the partition function 

Independent of Izergin and Korepin, in the case a = b (i.e. u = 0), an- 
other form of Z n was discovered conjecturally by Robbins in the context of 
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alternating-sign matrices (ASMs) and symmetry classes thereof, see [68]. As 
can be easily seen from Figures [5] and I ar weighted FPL configuration, 
where each straight loop segment is assigned a weight r, is equal to the gener- 
ating of weighted ASMs where each nonzero entry is assigned a weight r. Up 
to a simple factor, this is also the generating function A„(r 2 ) of r 2 -weighted 
ASMs of size n x n where each —1 is assigned a weight r 2 (each additional 
— 1 in an ASM also introduces an additional +1). The latter was conjectured 
by Robbins [68] to equal 

A 2n (r 2 ) = 2T„(r 2 )i? n _ 1 (r 2 ), ^ 2n+ i(r 2 ) = T n ( T 2 )R n {r 2 ). (13) 

where 

^ =j%(t(;-.^ (/ <») 

and 

(2n-l \ 

where 

= \2i + 1 - r) + \2i + + 1 +n - r) ' (16) 
The precise correspondence using the notation of the previous section is 

Z 2n = (sin 7 ) 2 "( 2 "- 1 ) A 2n (4cos 2 7 ), (17a) 
Z 2n+1 = 2 cos 7 (sin 7 ) 2 "( 2 " +1 )A 2 „ +1 (4cos 2 7). (17b) 

In fact, Robbins' conjecture was slightly more general and gave a generating 
function for refined ASMs. The generating functions R and T appear naturally 
in weighted enumerations of cyclically symmetric plane partitions [53]. 

The equivalence of the homogeneous limit of Izergin's determinant and 
Robbins' conjecture, i.e. equation (fT?]) , is only proved for r = 1 [42,80]. Ku- 
perberg and Robbins [43, 68] noticed several other such equivalences between 
homogeneous Izergin or Tsuchiy£0 type determinants and generating functions 
of the form (jT4j) or (p~5|) . Some of these were recently proved in [34] using a 
technique which seems immediately applicable to all the cases considered by 
Kuperberg and Robbins. 

3 Bond percolation, the 0(n = 1) model and the 
Razumov-Stroganov conjecture 

In this section we mention a (partially conjectural) relation between FPL 
diagrams and the 0(1) loop model. We will use this relation to generate FPL 

1 The Tsuchiya determinant is the generating function of horizontally or vertically 
symmetric FPL diagrams [76] 
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Fig. 11. Bond percolation clusters and O(l) cluster boundaries on a semi-infinite 
strip. Configurations are generated by repeated concatenation of double rows using 
the double-row transfer matrix. The particular boundary conditions chosen here are 
called closed or reflecting. 



statistics in a relatively easy way, without having to explicitly enumerate FPL 
diagrams. 

Imagine that each site of the square lattice is a reservoir of water. With 
probability p, water percolates between reservoirs along a bond of the square 
lattice. At p = 1/2, the model is critical, and equivalent to the dense 0(n = 1) 
loop model [12] on a square lattice. The loops of the 0(1) model describe the 
boundaries of the percolation clusters, see Figure [TT1 Many asymptotic prop- 
erties such as critical exponents of correlation functions can be computed for 
the 0(n) model using Coulomb gas techniques and conformal field theory, see 
Chapter 14 for an exhaustive overview. More recently, geometric properties of 
conformally invariant loops have been analysed using the stochastic Loewner 
evolution (SLE), see Chapter 15 of this book. 

Configurations of the 0(1) loop model can be generated using a transfer 
matrix, see Figure [TT] for the particular case of closed or reflecting boundary 
conditions. Schematically, the local blocks of the 0(1) transfer matrix are 
given by 

±v 2 + 
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The closed loops of the 0(1) loop model have weight n = 1. Loops ending on 
the boundary of the strip define a link pattern. For example, the link pattern 
corresponding to the bottom side of Figure [TT] has link pattern ()()()(). The 
transfer matrix T of the 0(1) loop model therefore acts on states indexed by 
a link pattern. 

3.1 The Razumov-Stroganov conjecture 

The largest eigenvalue of the transfer matrix of the 0(1) has eigenvalue 1. It 
was found in [5,32,59,64,65] that the corresponding groundstate eigenvector 
surprisingly is related to the statistics of FPL models. Denoting a link pattern 
by a and forming a vector space with basis elements \a), the groundstate 
eigenvector satisfies 



In the case of periodic boundary conditions, Razumov and Stroganov formu- 
lated the following important conjecture: 

The coefficient ip a equals the number of FPL diagrams with link pattern a. 



The RS conjecture generalises to other boundary conditions, in which case 
the eigenvector coefficient ip a of the corresponding transfer matrix enumerates 
symmetry classes of FPL diagrams, to be discussed below. This is explained 
in detail in [30]. The case that will be treated in most detail here is the 0(1) 
model on a strip, as in Figure 1 1 H for which ip a conjecturally enumerates 
horizontally symmetric FPL diagrams. 

Assuming the RS conjecture, we will use the 0(1) loop model to generate 
FPL statistics by solving (|19[) . and variants thereof for other boundary condi- 
tions. The particular boundary conditions we will use are periodic, cylindrical 
and closed. See e.g. [24,30,54,62,83] for other types of boundary conditions 
not considered here. 

Let us define the norm Ml of \*P) by 





(19) 



a 




(20) 



a 



and denote the largest element of \W) by ip max . The result of solving (|19[) for 
various boundary conditions is 
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Type 


Ml 




Periodic, L even 


A(L/2) 


A(L/2-l) 


Cylindrical, L even 


Awt{L) 


Aht{L — 1) 


Cylindrical, L odd 


A HT {L) 


A{(L-l)/2f 


Closed, L even 


A V (L + 1) 


Av{L) 


Closed, L odd 


A V (L + 1) 


Av(L) 



(21) 



where the numbers A, Aht and Ay are defined by, 

• The number ofnxn ASMs, 

4»>-nff$-w,42,... <m 

fc=0 v ' 

• The number ofnxn half turn symmetric ASMs, 

A HT (2n) = A(n) 2 JJ = 2,10,140,5544,... 

„-i = ° 2 (23) 

A„ T(2 „-l) = n !(<!§!!) = 1,3,25,588,... 

• The number of (2n — 1) x (2n — 1) horizontally (or vertically) symmetric 
ASMs, 

A v (2n-l)=n(3fc-l) g:^g:;j| = 1,1,3,2 6 ,64 6 ,.... (24) 
and its related version for even sizes (also denoted by N 8 in [13]), 

^v(2n)=n(3fc + l) ( i g|gg^ = 1,2,11,170,... (25) 

Mills et al. conjectured the number of ASMs to be A(n), which was proved 
more than a decade later by Zeilberger [80] and in an entirely different way 
by Kuperberg [42]. Kuperberg made essential use of the connection to the 
six- vertex model and its integrability. Conjectured enumerations of symmetry 
classes were given by Robbins [68] , many of which were subsequently proved 
by Kuperberg [43]. The properties and history of ASMs are reviewed in the 
book by Bressoud [13], as well as by Robbins [67] and Propp [61]. 
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3.2 Proofs and other developments 

The sum rules listed in Table [21] relating the norms ([20)1 of \<P} for differ- 
ent boundary conditions to symmetry classes of alternating-sign matrices, 
were originally obtained conjecturally. These sum rules have been proved alge- 
braically using an inhomogeneous extension of the transfer matrix, a method 
initiated and developed by Di Francesco and Zinn- Justin [19,21]. This has 
led to further interesting directions, not pursued here, such as the connec- 
tions between weighted FPL diagrams (or ASMs), plane partitions and the 
(/-deformed Knizhnik-Zamolodchikov equation [19, 20, 22, 23, 34, 57]. 

In an alternative interpretation, the 0(1) model is equivalent to a stochas- 
tic model defined on link patterns, the so called raise and peel model [33]. It is 
an open question how to define a stochastic model directly on FPL diagrams, 
by say the Wieland involutions G describe in Section 11.21 such that it has an 
equipartitc stationary state and reduces to the raise and peel model when the 
action of the operators O and U of Section 1X721 is divided out. Such a process 
would result in a direct proof of the Razumov-Stroganov conjecture. 

4 Symmetry classes of FPL diagrams 

We will now focus on FPL models defined on rectangular grids, corresponding 
to certain symmetry classes of square FPL diagrams. The two main reasons are 
that for such FPL models there is a natural boundary giving rise to additional 
structure, and that at the time of writing, for these models more results are 
known which are relevant to polygon models. 

4.1 Horizontally symmetric FPL diagrams 

For horizontally symmetric FPL diagrams (HSFPLs) one only has to consider 
the lower half of an FPL diagram. As explained in [30], due to geometric 
constraints one can further reduce the size of such half diagrams. Therefore, 
for L even, the reduced lower half of a horizontally symmetric FPL diagram 
of size (L + 1) x (L + 1) is an FPL diagram of size (L — 1) x L/2. The total 
number Zhsfpl(2?t.) of horizontally (or vertically) symmetric FPL diagrams of 
size (2n— 1) x n is known, and can be computed from the Tsuchiya determinant 
[43,76], 

Z HSF P L (2n) = A v (2n + 1) = f[(3k - ^ I (26) 

As can be seen from Table 1211 this number is equal to the norm Nm for 
the 0(1) model with closed boundary conditions and L — 2n. For odd system 
sizes, L — 2n + 1, the norm Ml equals the number of FPL diagrams of size 
L x (L — l)/2, which we will denote by Zhsfpl(2«. + 1). 
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L-l=7 

Fig. 12. Boundary conditions for an HSFPL diagram of size (2n — 1) x n = 7 x 4. 
The number of external terminals equals 2n = 8, hence the statistics of this diagram 
is generated from the 0(1) model with L — 8. 



innn 




£-1=13 

Fig. 13. An FPL diagram of size (L — 1) X L/2 = 13 x 7 with four nests. 



There are two interesting and natural statistics on HSFPLs which we will 
explain now. As noted above, to each FPL diagram is associated a link pattern. 
Each link pattern factorises in sets of completed links where, in terms of the 
parenthesis notation, the number of closing parentheses equals the number of 
opening parentheses. For example, 

(()())(((())())() = (()())•(((())())•()• 

Such completed links are called nests , and they provide a statistic for HSFPLs. 
An example of an HSFPL diagram of size 13 x 7 with four nests is given in 
Figure [H 

Another natural statistic is the number d* of loops connecting the leftmost 
loop terminals with the rightmost ones, i.e. loops connecting terminal i with 
2[L/2\ — i + 1 for i = 1, . . . , d*. It will be convenient to define d by 
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d = 



L - 1 



(27) 



where d is called the depth of an HSFPL diagram. An example of an HSFPL 
diagram of size (L — 1) X L/2 = 13 x 7 with three nests and depth d = 4 
(d* = 2) is given in Figure [Til 




d*=2 



L-l=13 



Fig. 14. An FPL diagram of size (L — 1) x L/2 = 13 x 7 with three nests and depth 
d = 4 (<f = 2). 



4.2 Depth-nest enumeration of HSFPLs 

In this section we will say that an FPL diagram is of size L, if it is of size 
(L — 1) x L/2 if L is even, or of size L x (L— l)/2 if L is odd. Let -P(X, d, m) be 
the number of such FPL diagrams of size L, depth d and having m + 1 nests. 
The nest generating function for diagrams of size L and depth d is defined by 

d 

V(L,d;z)=^2P{L,d,m)z m . (28) 

Let S(L, d) be the total number of HSFPL diagrams at a given size L and 
depth d. Obviously we have 

S(L,d) =V(L,d;l) = P(L,d + 1,0), (29) 

and 

^hsfplW = S(L, L VJ ) = HL, L^J ; I)- (30) 

Based on the RS conjecture, Mitra et al. and Pyatov have conjectured the 
exact form of S(L,d) [54,62]. Here we give this conjecture in the following 
form: 
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Conjecture 1 The total number of HSFPL diagrams at a given size L and 
depth d is given by 

j-j 2 fc (i/2) fc r(L-2fc + i) rC L ~ 6 k+3 )r( 2L - 6 k + 6 ) ' 1 ' 

Assuming the RS conjecture, the formula for S(L, d) has recently been proved 
[34]. 

Pyatov also found an exact formula for P(L,d,m) [62] which fits exact 
data for small system sizes (L < 18). He conjectured that this formula holds 
for all L, d and m. In terms of the nest generating function this conjecture 
can be stated as follows. 

Conjecture 2 The nest generating function is given by 

V(L, d: Z ) = S(L, d - 1) 3 F 2 ( ~ d ll - 2 f L - d d + 1 * 5 4*) . (32) 

Note that Conjecture Q] follows from Conjecture [5] due to the evaluation 

(-d,L-2d,L-d+\ A _ S(L,d) 
3F \ -2d,2L-2d+l ' 4 J- S{L,d-lY (33) 

which is a consequence of one of the strange evaluations of Gessel and Stanton 
[28]. For d = [(L — 1)/2J , the formulas in Conjecture Q] and Conjecture [2] were 
given in [30]. 

By convention, P(L,d,m) have the following boundary values: 

P(L, d, m = -1) = P(L, d,m= -2) = P(L, d, m = d+ 1) = , (34) 

and we also note the boundary condition 

P(L,d,m = 1) = (L-2d)S(L,d- 1) . (35) 

It was found in [1] that the function P(L,d, m) is completely determined by 
these boundary conditions and the following interesting bilinear relation called 
the split hexagon relation, 

P(L + l,d+ l,m)S(L - l,d - 1) 

(36) 

= P(L - 1, d, m)S{L + 1, d) + P{L, d-l,m- 2)S(L, d+l). 

Summing up over m = 0, 1, . . . d + 1 in (|36p reproduces the hexagon relation, 
or discrete Boussinesq equation, for S(L,d), see [62]. 
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Cyclically symmetric transpose complement plane partitions 

Somewhat outside the scope of this book, we note the following interesting fact 
observed in [34]. The total number of nests at a given depth, S(L, d), is equal 
to the number of punctured cyclically symmetric transpose complement plane 
partitions [14], see Figure fT5l This can be seen by enumerating the number 
of non-intersecting lattice paths in the South-East fundamental domain of 
the plane partition. Using the Gessel-Viennot-Lindstrom method [29,49] one 
obtains a determinant of the type (|14[) with r = 1, which can be evaluated in 
factorised form [14]. This form equals the expression in (13 1 1) . 




4.3 Average number of nests in HSFPL diagrams 

The average number of nests in HSFPL diagrams at depth d and size L, 
denoted by, (1 + m)d*, is defined as 
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1 



d 



(l + m) d . 



Zhsfpl(L) 



^T(l + m)P(L,d,m) 




(37) 



For notational clarity we will suppress the dependence of (1 + m)d* on L and 
recall that 



With the data P(L,d,m) we can calculate (m) d ,: 



The Mathematica implementation of the Gosper-Zeilberger algorithm 
[35,78,79] by Paule and Schorn [58], is able to recognise (m)d* in an almost 
factorised form. Let L = 2n, then define ji n {d*) by 



The expression /i n (d*) turns out to be summable in factorised form, giving 
rise to 



Ms- =-fo-2d)+2 2 /3 ; pL u^; P ^; f ,y - ( 4o ) 



This formula also holds for odd values of L. 
4.4 Half-turn symmetric FPL diagrams 

In the case of half turn symmetric FPL diagrams (HTSFPLs) it also suffices to 
consider only the lower half of an FPL diagram, but the boundary conditions 
on the top row of the half diagram are different from HSFPLs, see Figure [THJ 
The total number of HTSFPL diagrams is given by [43] 




d 



logV(L,d;z). 



(38) 





Z H T(2n) = A H T(2n) = 2]J 



3(3fc + 2)!(3fc- l)!fc!(fc- 1)! 
4(2fc + l)! 2 (2fc- l)! 2 



(41) 
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L-l=7 

Fig. 16. Boundary conditions for an HTSFPL diagram of size (2n — 1) x n = 7 x 4. 
The arcs at the top are additional edges which may contain loop segments. The 
number of external terminals equals 2n = 8, hence the statistics of this diagram is 
generated from the periodic O(l) model with L = 8. 



Care has to be taken when defining link patterns and nests for HTSFPL 
diagrams. External terminals can be connected in two distinct ways depending 
on whether the corresponding loop runs over an odd or even number of the arcs 
on the top of the diagram. In the case of an odd number of arcs, we exchange 
the parentheses denoting the connection of a pair of sites. For example, the 
connectivity of the HTSFPL diagram in Figure [17] is denoted by 

)•()•()•((), 

where the dots again denote the factorisation of link pattern into nests. Fig- 
ure [T7] thus denotes an HTSFPL diagram with three nests. 

As in the case of horizontal symmetry, there exists a conjecture for the 
nest distribution function [30], but in this case only for L = 2n and d* = 0. 
Let P(L, m) denote the number of half-turn symmetric FPL diagrams with 
77i + l nests, and define the nest generating function by 

n-l 

]TP(L,m)z™ ( 42 ) 

m=0 

Conjecture 3 The nest generating function for half-turn symmetric FPL di- 
agrams is given by 

P(2n;z) = Z H T SF PL(2n) i - J - T 3^( v 2 _ 2n 2 + 2n ;4z). 

The average number of nests in HTSFPL diagrams of size L = 2n having 
1 + m nests, denoted by (1 + m), is defined as 
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L-l=7 

Fig. 17. A HTSFPL diagram of size (2n- 1) x n = 7 x 4 with link pattern )•()■()•((), 
having three nests. 



n-1 



(1 + m) 



V (l + m)P(L,d,m) = 1 + z— logV(L;z) (43) 



•^HTSFPL (L) ^ Q 



Knowing the nest generating function we may compute (1 + to), which turns 
out to be summable [30]. 

Conjecture 4 The average number of nests in HTSFPL diagrams of size L, 
is given by 



(1 + to) = n J \ 



3j + l 
L : 3j + 2' 



5 Phase transitions 

5.1 Bulk asymptotics and phase diagram 

The phase diagram of the FPL model can be derived from the asymptotics of 
the partition function Z n defined in (|10p . The leading asymptotics of Z n for 
general values of r has been computed by Korepin and Zinn-Justin [41] using 
the Toda equation satisfied by a n [74] , 

d 2 cr„ /dcr n \ 2 

an ~j — 5~ ~ ~3 — = <Tn+i<r n -\- (44) 
au z \ d u J 

Writing (7 fi cLS Bi matrix model integral [82], further subleading asymptotics 
were computed by Bleher and Fokin [10] and Bleher and Liechty [11] using 
orthogonal polynomials. For special values of 7 this method was first employed 
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by Colomo and Pronko [17]. The final result for < t 2 = 4 cos 2 7 < 4 (1 > 
A > —1) is that for some e > 0, 

Z n = Cn K exp[fn 2 (l + 0(n- 6 ))] , (45) 

where C is a constant and 

7r sin(7 + u) sin(7 — u) 



27 cos(7ru/27) 
1_ 2 7 2 
12 3vr(7r - 27) ' 



(46) 
(47) 



The result (|45|) is valid in the so-called disordered (D) phase < r 2 < 4. 
There is a phase transition to an ordered phase at t 2 = 4 where the vertices 
with weight c are favoured and the perimeters of the polygons in the FPL 
model consist of elongated straight lines. In terms of the six-vertex model 
this is called the anti-ferromagnetic (AF) phase. At r = 0, i.e. a — b + c or 
b = a + c, there is another phase transition to a so-called ferromagnetic phase, 
where, respectively, the a- or 6-type vertices dominate. The complete phase 
diagram is given in Figure [TBI 



b/c 




1 a/c 



Fig. 18. Bulk phase diagram of the FPL model. The phases are traditionally 
called disordered (D), ferro-electric (F) and anti-ferro-electric (AF), cf. the six- vertex 
model. The arc corresponds to the free-fermion condition A — (t 2 = 2), and the 
line a = b corresponds to r 2 -enumerations of ASMs. On this line r 2 = c 2 /a 2 = 
2(1 — A), and the D-AF phase transition takes place at r 2 = 4. 



The phase diagram and the Bethe-Ansatz solution of the six- vertex model 
for periodic and anti-periodic boundary conditions are thoroughly discussed 
in the works of Lieb [45-47], Lieb and Wu [48], Sutherland [73], Baxter [7], 
and Batchelor et al. [3]. 
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5.2 Asymptotics for symmetry classes at r = 1 

In this section we determine the asymptotics of equally weighted horizontally 
and half-turn symmetric FPL diagrams for r = 1 or 7 = 27r/3, correspond- 
ing to the numbers given in Table 1211 The leading asymptotic form of these 
numbers, which are all products over factorials, can be computed using the 
Euler-Maclaurin approximation. Full asymptotics can easily be derived using 
Barnes' G-function [2], which satisfies 

G(z + l) = r(z)G(z), G(l) = l, (48) 

and whose leading asymptotic behaviour is given by (see e.g. [56]), 

log(G(z + 1)) = z 2 Q log z - ^ + \z log 2tt - 1 log z + 0(1). (49) 

In the case of A(n) and Aht(^), a detailed asymptotic analysis including the 
lower order terms was carried out by Mitra and Nienhuis [55] . Here we list only 
the leading asymptotics of the FPL numbers relevant to the current context. 
The generic asymptotic form of the numbers is 

log Z L = s Area + f Surface + x log(Length) + 0(1), (50) 

where the bulk and boundary entropies are given by 

The critical exponent x is a universal quantity. In detail, the cases relevant 
for this chapter are 

• FPL diagrams, L even 

The number A(L/2) counts FPL configurations on an L/2 x 1/2 square 
grid, the area of which is \L 2 ■ We thus find, 

logZ(L) = log A(L/2) = ^s L 2 - A \ ogL + (1). (52) 

• Half turn symmetric FPL diagrams, L even 

Awt{L) counts the number of FPL configurations on half an L x L square 
grid, the area of which is \L 2 . We thus find for L even, 

logZ H T(£) =logA HT (i) = \s L 2 + 4 log L + 0(1). (53) 

Z lo 

• Half turn symmetric FPL diagrams, L odd 

Awt{L) counts the number of FPL configurations on a square grid of 
dimension L x (L — l)/2, the area of which is \L{L — 1). We thus find for 
L odd, 

logZ HT (i) = log Aht(L) = ls L(L - 1) + ^ logL 2 + 0(1). (54) 

2 6b 
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• Horizontally symmetric FPL diagrams, L even 

Ay(L + 1) counts the number of FPL configurations on an (L — f ) x L/2 
rectangular grid. We find, 

logZ HS FPL(i) =log^ v (£ + l) - ^s L(L-l) + f L-^logL + O(l). 

(55) 

• Horizontally symmetric FPL diagrams, L odd 

For L odd, Ay(L + 1) counts the number of FPL configurations on a 
L x (L — l)/2 rectangular grid. We find, 

logZ HS FPL(i) =logA V (£ + l) = i« i(i ~ 1) + /oi + ^logi + 0(l). 

(56) 

Note that because the upper boundary for FPL diagrams corresponding to 
HSFPLs is not fixed, see e.g. Figure [T2l there is a nonzero boundary entropy 
in logZ H sFPL(i)- 

5.3 Nest phase transitions 



From Section 14.31 we recall that the average number of nests is given by 
+ where 

(m)% = z^-logV(L,d;z), (57) 
az 

with V(L,d;z) given in Conjecture [21 The asymptotics for (m)^» as L — > 
oo can be derived from the hypergeometric equation satisfied by V(L,d]z). 
Taking L = 2n this gives 

9(9 + 1 + 2d* -2n)(9 + 2 + 2d* + 2n)P(2n, d; z) = 

Az(9 + 2 + d*)(9 + I + d* - n)(6 + 3/2 + n)V(2n, d; z), (58) 

where 9 = zd/dz and 

L—l 
2 

Assuming that d* —0(1), we discriminate the cases z < 1, z = 1 and z > 1 . 
• z < 1 

In this case, up to an overall constant factor, the leading asymptotics of 
V(2n, d; z) will be polynomial in n. Neglecting lower order terms, equation 



d*=n-l-d*. (59) 



([58)1 reduces to, 

9V(2n, d; z) = z{9 + 2 + d*)V(2n, d; z). (60) 
We thus we find (1 - z)V'(2n, d; z) = (2 + d*)V(2n, d; z) and 

(m)3. = (2 + d*) T - L - (»-°o). (61) 
1 — z 
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• z = l 

In PO]) an exact expression was given for (m) d , at z = 1. Asymptotically 
we find that for L - 2d = 0(1), 



L-2d+l ' 



{™n-™-j^=&fL*/* + 0(l), (62) 
which for L = 2n can be written as 

7-./ 2d*+3 \ 

Hf ^jyj t^ + Otl), (63) 

• z > 1 

In this case, and when d is of order n, the leading asymptotics of "P(2n, <i; z) 
will be of the form p(n)z n , where p(n) is a polynomial in n. This means 
that QT '(2n, d; z) is of the same order as n'P(2n, d; z) and (|58[) reduces in 
leading order to 

(6» 3 - 4n 2 6»)P(2n, d; z) = 4z(6» 3 - n 2 9)V(2n, d; z). (64) 

Using (|57p one can derive the following equation for {m) c d , , 

4n 2 (z - l)(m) c d , = (4z - 1) (9 2 (m) d , + 3(m)%0(m)% + ((m) d ,) 3 ) , (65) 

which in leading order when (m) d , ~ n reduces to (4z — 1)((™)^*) 2 = 
4n 2 (z - 1) and thus 



|L_J_L + 0(1). (66) 



The scaling behaviour near the phase transition at z = 1 is governed 
by a single exponent, the cross-over exponent <f> [9], On general grounds one 
expects, 

1)V*-1 ( z > 1) 

(z = l). (67) 
(1-z)- 1 (z<l) 

Indeed, we find such scaling behaviour for {m) d * with (f> — 2/3. 



Scaling function 

In [33] an analysis has been carried out to obtain the nest scaling function for 
L = 2n and d = n — 1, i.e. d* = 0. Following Polyakov [60], we expect the 
following scaling form of the nest distribution function at the critical point, 
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P(2n,n-l,m) 1 / 1 + m , 

j 7i — I r (n->oo), (68) 



S(2n,n-1) (l + m)o \(l + m) 

where (1 + m) = 1 + (fn)o- The large x behaviour of f(x) is related to the 
exponent <fi [15], 

lim /(a) ~x s e~ ax \ S = — !— , (69) 

rr^oo 1 — 

where a and s are constants. The behaviour of /(a;) for small x is related to 
the large n behaviour of the probability P(2n, n — 1, m)/S(2n, n — 1), 

lim /(*) = &«:* =* b = lim lim (1 + (rn) ) 1+ " P g"' " - ^ , (70) 



from which we find 



tf = 1 ' 6= TW (71) 



Assuming that the full scaling function is of the form x" 9 e axS , for all values 
of x, and using the normalisation condition 

f(x)dx = l, (72) 

we find that 

f(x) = bx e- bx3/3 . (73) 

In Figure [TH] we compare the scaling function (|73p with a numerical evaluation 
of ([68]) for n — 300. When 8* > 0, the value of (f>, and hence that of S, is not 
changed, but it follows from ([70]) that the value of the exponent $ changes to 

tf = l + 2d*. (74) 

The full scaling function is not known in this case. 



5.4 Half turn symmetry 

The following analysis closely follows that of the previous section. We are 
interested in the asymptotics as n — > oo of the average number of nests defined 
in (|57| . This can be inferred from the hypergeometric equation for V(2n\ z), 

6{6+l+2n)(6+l-2n)V(2n; z) = 4z(6+3/2)(8+l-n)(9+l+n)V(2n; z). (75) 

Again we discriminate the cases z < 1 , z = 1 and z > 1 and remind the reader 
that L = 2n. 



Fully packed loop models on finite geometries 



25 



fix) 1.0 




0.0 0.5 1.0 1.5 2.0 

Fig. 19. The scaling function f(x) defined in (|73[) compared to a numerical evalu- 
ation (dots) of (J6H1) for L = 2n — 600. 

• z < 1 

In this case, up to an overall constant prefactor, V(2n;z) will grow as a 
polynomial in n and, neglecting lower order terms, (|75p reduces to, 

6V(2n;z) = z(6 + 3/2)T(2n; z), (76) 

so that we find (1 — z)V'(2n; z) = |7 , (2n; z) and thus 

(l + m)«^ + 0(l). (77) 

• z = l 

For this case, and exact expression was given for (m + 1) in Conjecture |H 
Asymptotically we find 

(1 + m) = n if ^ « ^ ^ 2/3 + <W (78) 

• z > 1 

Here, up to an overall constant prefactor, 7 , (2n;z) will grow as p(n)z n 
where p(n) is a polynomial in n. This means that 9V(2n; z) will be of the 
same order as riP(2n\ z) and (|75[) reduces in leading order to 

(6> 3 - 4n 2 6»)P(2n; z) = 4z(6> 3 - n 2 (9)7>(2n; z), (79) 

which is the same as We thus find that 

< 1+ ™> + ( 8 °) 




26 Jan de Gier 



For half-turn symmetric FPL diagrams we find the same cross-over exponent 
<j> = 2/3 as for horizontally symmetric FPL diagrams. 

6 Conclusion 

We have described a model of tightly packed, nested polygons on the square 
lattice. We hope that the study of such tightly packed polygons is relevant 
to other polygon models described in this book. The advantage of the model 
described in this chapter is that many exact results can be obtained, even on 
finite geometries, due to its relation with the exactly solvable six-vertex and 
0(n = 1) lattice models. In particular, the statistical mechanical partition 
function can be obtained rigorously on finite square patches of the square 
lattice. The free energy can then be obtained analytically and in the ther- 
modynamic limit. The fully packed loop model undergoes a well known bulk 
order-disorder phase transition as a function of an anisotropy parameter as- 
sociated to the straight segments of the polygon boundary. 

We have furthermore shown that it is possible to obtain closed form expres- 
sions for partition functions of important subsets of fully packed loop configu- 
rations. Two examples of such subsets are horizontally symmetric fully packed 
loop models of depth d, and half-turn symmetric fully packed loop models of 
depth d. These closed form expressions have been obtained experimentally, 
and remain conjectures at the time of this writing. In addition, using hyper- 
geometric summation identities, we were able to compute the average number 
of polygon nests at the boundary in closed form. Asymptotic analyses allowed 
us to study a boundary phase transition as a function of the nest fugacity, 
and we obtained a crossover exponent <fi — 2/3. At criticality, we derive a 
scaling form for the nest distribution function which displays an unusual non- 
Gaussian cubic exponential behaviour. 

To conclude it should be said that while some of the exact results presented 
in this chapter are more than what one would hope for from a physicist's per- 
spective, where numerical techniques are often all what is available, they are 
just the starting point for a mathematician. Although great progress has been 
made in recent years in understanding fully packed loop models, proving con- 
jectures such as the nest distribution function and many related combinatorial 
results, remains a fascinating and completely open problem. It is for reasons 
such as these that polygon models in all shapes and sizes will continue to 
inspire future research. 
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